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MOTIONS OF AN AIRCRAFT IN RESPONSE TO STEP CONTROL 
DEFLECTIONS AND TO THE INFLUENCE OF 
PROPORTIONAL AUTOMATIC CONTROL 
By Howard J. Curfman, Jr. 


SUMMARY 


A study has been made of the effects of two nonlinear stability 
derivatives, caused by the nonlinear variations of pitching-moment 
lift coefficients with angle of attack, on the longitudinal motions of 
an aircraft. Theoretical methods involving the Laplace transformation 
and the procedures of nonlinear mechanics have been presented along with 
electrical-analog results for the responses of a canard aircraft to 
step control deflections and to the influence of two types of propor- 
tional automatic control. 

In all cases except the example illustrating the. procedures of 
nonlinear mechanics, the nonlinear variations have been approximated by 
three linear segments i The principal means of studying the effects of 
the nonlinearities were through the characteristics of the time responses 
in angle of attack and, in some cases, pitching velocity. The changes 
in the period and damping of the oscillations due to the nonlinearities 
have been discussed. The effects of the autopilot proportionality con- 
stant on system stabilit^y were also investigated. The occurrence of 
continuous hunting oscillations was predicted and demonstrated for the 
attitude stabilization system with proportional control for certain 
nonlinear pitching -moment variations and autopilot adjustments. 


INTRODUCTION 


In the classical treatments of the study of aircraft dynamics an 
assumption, is usually made that the forces and moments resulting from 
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small disturbajices in position, velocity, and acceleration are linear 
variations with these quantities. However, some aircraft configurations 
do not exhibit linear force, and moment characteristics in some flight 
conditions and, therefore, the linearized methods cannot be a satis- 
factory approximation over an appreciable range of the parameter 
involved. This paper is concerned with the problem of aircraft dynamics 
as affected by nonlinear stability derivatives and, in particular, with 
the short -period longitudinal mode of motion of an aircraft flyirfg at 
constant velocity. Consideration is also given to the problem of 
applying simple forms of automatic control to such an aircraft; 

The present work is presented in two independent parts. The first 
part considers several analytical approaches that may be used in the 
analysis and the study of the transient motions of the aircraft. These 
methods include the use of the operational calculus (Laplace transforma- 
tion), the stability criterion of Routh, and the theories of nonlinear 
mechanics. The second part presents some of the results obtained by 
using electrical analog equipment. The cases presented are for a 
hypothetical canard aircraft assumed to be flying at- supersonic speed. 


SYMBOLS 


a angle of attack, positive when nose is above relative wind 

vector, radians unless otherwise noted 


d 


da 

dt 


aircraft attitude or pitch angle, positive when nose is above 
horizontal, radians 


A _ de 
® - dt 


6 aircraft elevator deflection angle, positive when trailing 

edge is down, radians unless otherwise noted 

V aircraft velocity, feet per second 

q dynamic pressure, pounds per square foot 

S wing (reference) area, square feet 

c mean aerodynamic chord, feet 

mass of aircraft, slugs 


m 
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ly pitching moment of inertia, slug-feet square 

D differential operator (d/dt) 

t time 

Cm pitching -moment coefficient (Pitching moment/qSc) 

Cl lift coefficient (Lift/qS) 


°2V 


" 1 %. 




m 

0 £ 

2V 


Cm5 




Ct - 

- 56 “ 

kp slope of a linear segment representing a portion of nonlinear 

Cm(o,) function, per radian 

^2 slope of a linear segment representing a portion of nonlinear 

, - CL(a) function, per radian 

a,b,c,d constants used in equations (6) and (9) and defined in appendix 

o,Ti constants used in equations (2) and (3), respectively 

undamped natural frequency of aircraft alone, radians per 
second 

C nondimenslonal damping ratio of aircraft mode of motion 

s Laplace transform variable 

K general autopilot constant, ratio of control surface deflec- 

tion to applied error angle 
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If - 5 

^2 - 0 

tti reference or desired angle of attack for closed-loop system 

with angle -of -attack feedback, radians 

9-1 reference or desired attitude angle for closed-loop attitude - 

stabilization system, radians 

constants required in cubic representation of nonlinear Cui(a) 
function (equation (12)) 

slope of trajectory in phase plane example (d0/da) ' 


ANALYSIS 


The analysis of the effects of nonlinear stability derivatives on 
the longitudinal motions of a canard aircraft is presented in three sec- 
tions: the first is a statement of the problem and the assumptions 
involved; the second is a theoretical or analytical solution for the 
required, transients and the study of stability; and the third is a dis- 
cussion of the methods of nonlinear mechanics for these dynamical 
studies. 


Statement of Problem 

The purpose of this paper is to present methods for determining the 
nature of the transient responses of aircraft having nonlinear stability 
derivatives and to show some qualitative results of preliminary 
investigations. 

In the present paper, consideration is given to the longitudinal 
motions of a canard configuration with the angle -of -attack transients 
resulting from step elevator deflections serving to illustrate the 
effects of the nonlinear stability derivatives. Two' longitudinal - 
stability derivatives were considered to have nonlinear characteristics 
which were caused by the variations of pitching-moment coefficient and ' 
lift coefficient with angle of attack. The main emphasis herein was 
placed on the effects of the nonlinear variation of pitching-moment 
coefficient with angle of attack since the slope of this curve is a 
measure of the static stability of the aircraft. Figure 1 illustrates 
the nature of the nonlinearities under discussion. Such nonlinearities 
would be apparent from wind-tunnel tests. 


X 
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In the analysis and study of the aircraft longitudinal motions, 
two degrees of freedom with constant forward velocity and disturbances 
from level flight have been assumed. This assumption results in a con- 
sideration of the short-period mode of motion only; therefore, only the 
effects of the nonlinearities on this mode will be apparent. 

In addition to studying transients of the aircraft in response to 
step elevator inputs, consideration was also given to the response 
characteristics resulting when such an aircraft is subjected to simple 
forms of automatic stabilization. Two cases were investigated in this 
respect and are illustrated by block diagrams in figure 2. In both 
cases the automatic control consists of an error-sensing device and a 
zero-lag proportional servomotor. In figure 2(a) angle-of -attack sta- 
bilization is considered with the control-surface deflection S ^being 
proportional to the difference between a reference angle of attack ai 
and the actual angle of attack a. In figure 2(b) attitude stabiliza- 
tion is considered where the control -surface ‘deflection .6 is propor- 
tional to the difference between a reference attitude angle G± and 
the actual attitude 0. 


Theoretical Methods 

Equations of motion .- The longitudinal equations of motion of an 
aircraft having two degrees of freedom with constant forward velocity 
and disturbances from level flight in terms of the differential oper- 
ator D = -^ are 
dt 


qSc ® ® 2V ~ Cm(a) - - Cjng6 


g D0 . § Da . Ci,(a) - Cl^5 


( 1 ) 


qS qS 


The terms which are nonlinear functions of angle of attack are written 
as Cm(a) and CL(a) and represent, respectively, the pitching-moment 
coefficient and lift coefficient that exist at the instantaneous angle 
of attack. 


Based on the usual assumptions , of linear variations of Cm nni^ Cj, 
with a, equations (l) remain as ordinary liriear differential equations 
with constant coefficients and are handled by familiar mathematical 
methods. The types of nonlinearities illustrated in figure 1, however, 
obviously make linearization of these parameters over the entire range 
of angle of attack a poor approximation. Since the methods for the 
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solution of linear differential equations are so well known, an approxi- 
mation of the nonlinear functions with linear segments, the equations of 
which can be written for, a specified range of angle of attack, seems 
advisable. Figure 3 is an illustration of this method of representing 
a nonlinear function. The three straight lines shown are used to repre- 
sent a variation of pitching-moment coefficient with angle of attack 
similar to the one shown in figure 1, where the slope at outer a is 
twice the slope at central a. Each segment of this representation 
canbe expressed as follows; 




where ki arui a must be specified for each range of angle of attack. 
In general, therefore, the representation of the nonlinearities are 
written as 


Cj„(a) = kj^a + a (2) 

CL(a) = k2tt + Ti (3) 

where the values of k^, a, and r\ are constants and apply only 

for a specified range of angle of attack. 


The substitution of equations (2) and (3) into equations (1) 
results in linear differential equations with constant coefficients, 
namely. 


^ ^ D0 - kia - ^ Da = CmgS + a 


g D0 - ^ Da - k 2 a = q 




( 4 ) 


These equations apply only to the region of angle of attack in which the 
given constants are as defined. 


Aircraft transients .- Since variations in angle of attack following 
step elevator deflections may be in regions where the definition of the 
nonlinear functions has changed, complete knowledge of the motion is 
required and solutions of the equations must be written in terms of 
necessary conditions existing at an arbitrary time. In other words, 
the values of the angles and their rates of change must be known at the 
time one set of equations ceases to apply so that the same conditions 
and, therefore, continuity may be satisfied when the next- set of 
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equations becomes applicable. ’With these considerations the use of 
operational methods based on the Laplace transformation affords a means 
of handling the difficulties already noted and provides the required 
analytical expressions for the desired transient responses. Refer- ' 
ences 1 to 4 give adequate discussions of the use of the Laplace trans- 
formation in the solution of differential equations. Since the usual 
forms of the Laplace transformation require the initial conditions at 
zero time, it is convenient to make solutions for the required motions for 
each range of angle of attack and it may be necessary to redefine the time 
as zero whenever the equations must be changed because of the repre- 
sentation of the nonlineaxities. Such a procedure is obviously time 
consximing and tedious since solutions for all the parameters are 
required for a complete definition of initial conditions. A rather com- 
plete discussion of this procedure can be found in reference 5. 

Closed-loop stability .- The stability of a complete closed-loop 
system consisting of the aircraft and automatic -control unit, such as 
illustrated in figure 2, has been theorized by many authors when the 
system is linear throughout; see, for example, references 2, 6, 7, and 8. 
Briefly, this procedure is to plot the frequency response (reference 9) 
of the aircraft as a Nyquist diagram (references 2, 7, and 8) or as an 
amplitude -phase plot (references 2 and 6) and then to determine the • 
autopilot constant K required to make the system stable and give, as 
closely as possible, the desired response. _ Since the frequency-response 
method requires a completely linear system, such a method cannot be 
applied in the conventional manner and have meaning because of the 
definition of the nonlinear functions. 

\ 

This stability, however, can be examined by the methods of Routh 
and Hurwitz described in references 10, 11, and 12. This proced\ore 
requires the characteristic equations that exist for each region used 
in the representation of the nonlinear functions. If these investiga- 
tions indicate stability in all regions of the parametric representa- 
tion, the stability of the closed loop is generally assured. Such 
methods could also be used, of course, to determine the values of the 
autopilot constant K required to assure con5)lete stability. 

The application of this technique for the cases of angle-of -attack 
feedback and attitude stabilization is as follows: 

(a) Angle -of -attack feedback: The control equation is written as 

follows for the angle-of -attack feedback system (see fig. 2 (a)): 

6 = Ki(ai - a) (5) 

When equation (5) is substituted into equations (!<■) and the Laplace 
transformations of the resulting expressions are formed, the solution 
for a can 'be written in terms of the transform variable as 
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a(s) = 


K3_(as^ + bs)a^(s) + |^.C.(s^ 

! + ( 2^ (Du + Kxa ) s + + Kib 


( 6 ) 


where the term |T.C.(s^ is the function of s, involving the 

initial, conditions and constants, that results when the control equation 
and equations (4) have been transformed and solved for a(s), and the 
other coefficients axe constents. These expressions are tabulated in 
terms of the aircraft peurameters in the appendix. 


The stability criteria of the characteristic equation of (6) are 

2^0ijj + Kja > o'] 


and 


CD 


+ K^b > 0 


J 


(7) 


The inequalities (7) allow a thorough investigation of the effects of 
Ki on stability. These criteria must be checked qf course in each 
range of a used in defining the nonlinear functions which appear in 
equations (1) if the type of- motion is to be completely investigated. 
The types of motion possible are similar to those of a simple spring, 
mass, viscous damping system which also has a quadratic characteristic 
equation (see, for exaaqjle, references 2 and 13). The equation which 
describes the attitude variations as this type of system operates can 
be determined by . solving for 0(s) after equations (i)-) and (5) are 
transformed as described previously. 


The angle -of -attack feedback system does not exhibit a zero error 
at steady state even for the completely linear case. This condition 
may be seen by applying the final-value theorem (see, for example, 
reference 2, pp. 72-73) to equation (6) and noting that the steady- 
state angle of attack does not equal the input and, hence, the 

error - a does not tend to zero. This characteristic also is 
present when the nonlinear functions are considered'. For the nonlin- 
earities considered herein, all a transients are completely stable, 
that is, exhibit steady- state values. 


(b) Attitude stabilization: The control equation is written as 

follows when attitude stabilization is considered (see fig. 2(b)): 

6 = K2(0i - 6) 


( 8 ) 
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When equation (8) is substituted into equations ( 4 ) and the transforma- 
tions are made, the complete closed-loop solution for 0(s) becomes. 


K2(cs + d)0i(s) + [l.C. (sH 
~ — * 0 

s3 + 2^01^3^ + + K2 c)s rf K2d 


( 9 ) 


where 


(i.c.(sri^ 


is the function of 


involving initial conditions 


and constants, ’that results when the control equation and equations (V) 
have been transformed and solved for ^(s), and the other factors are 
constants as defined in the appendix. 


The criteria for stability of the cubic characteristic equation of 
equation (9) require that all coefficients be positive and that 


(2Co)n)(cUn^ + K2 c) - K2d > 0, (10) 


These criteria allow a complete investigation of the effect of on 

stability and should be checked for each angle -of -attack range where 
definition of the nonlinear functions changes. The angle -of -attack 
variation a(s) for this system can be solved for from the transformed 
equations ( 4 ) after the -control equation (8) is considered. 

For the completely linear case when stability is assured, equa- 
tion (9) reveals that the system is a zero -error system, the output 
tending to equal the input and the error ( 0^ - 0) tending to . equal zero. 
For the nonlinearities considered herein, two conditions affecting 
operation of the system can exist, and these conditions result primarily 
from the nonlinear Cj„(a) function. (For normal C^Ca) functions where . 

the slope remains positive and is not greatly changed in the angle -of - 
attack regions, this nonlinearity will not affect the stability or 
oscillatory characteristics to any appreciable extent.) For the system 
having a nonlinear Cj^(a) function and a satisfying the stability 

criteria for all values of a, it may be shown from the examination of 
steady-state a and steady-state 0 following a step input ■ 9± that 
the system remains a zero-error system. Consider secondly the case 
where the CmCa) function has an unstable slope at central a and a 
stable slope at outer a. and a K2 such that stability is assured at 
outer a , but not at central a. The instability at central a means 
that the angle -of -attack motion in this region is divergent to outer a; 
that is, no steady-state a. can exist at central a. Upon examination 
of a(s) it will also be noted that no steady-state a can exist at 
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outer a; however, the motion at outer a is stable, and the a tends 
to the region of central a where the motion has been shown to be 
divergent. This process is summarized as follows: 

(1) At central a, the angle-of -attack motion is always divergent 
to outer a 

(2) No steady-state a exists at outer a 

(3) The tendency of the motion at outer a is toward the region 
of central a 

It is concluded, therefore, that the angle -of -attack motion should, 
eifter some initial transients, become some form of continuous hvinting 
oscillation (nonsihusoidal) rendering the complete system continuously 
oscillatory. A constant -amplitude constant -period motion would also be 
suggested from the summary of the angle -of -attack tendencies when exam- 
ined on an energy basis, an energy balance becoming established between 
that energy fed in by the unstable mode and that dissipated in the 
stable regions. 

In summary, the attitude-stabilization analysis reveals the fol- 
lowing significant facts: 

f 

(1) The autopilot proportionality constant K2 can be adjusted 
for complete stability at all values of a and renders the system a 

zero-position-error system for the nonlinearities considered herein. 

/■ 

(2) If the CiB(a) nonlinearity is such that the constant K2 
assures stability at outer a but not at central a, continuous 
hunting oscillations can exist. 

Closed-loop transients .- Equations (6) and ( 9 ) are completely 
expressed when ai(s) and ^9i(s) are known. For a step, function of unit 
magnitude. 


0j^(s) = aj,(s) = j 

Upon substitution into the proper equations, the transform fxmction is 
known and the inverse transformation to a(t) and/or 0(t) can be 
determined . 


Methods of Nonlinear Mechanics - 

For the types of nonlinear functions represented in figure 1, the 
use of a cubic equation is suggested as a possible representation of the 



NACA RM L50L11. 


11 


nonlinearity over an adequate range of angle of attack. Such a cubic 
representation would, of course, make the differential equations of 
motion of the aircraft nonlinear. To many such problems, however, the 
methods of nonlinear mechanics have proved a powerful means of satis- 
factory analysis, although the results are seldom available as the 
familiar time responses. References Ik and 15 afford very complete 
discussions of the methods and techniques used in applications of the 
theory of nonlinear mechanics. 

For the present paper a single application of the methods of non- 
linear mechanics will be discussed in detail to illustrate the applica- 
tion and usefulness of such methods to the problem considered. For the 
case chosen and, in general, for the cases where the principles of non- 
lineax mechanics are readily applicable, the proce;dure is to rewrite 
the differential equations of motion into an independent system, that 
is, into a system of equations in which time appears only as a differ- 
ential dt. Such a representation is possible for many cases or can be 
derived by a suitable change in variable. After such a revision, a 
phase plane can be defined and the characteristics of the motion can be 
determined by trajectories in this phase plane. The coordinates of the 
phase plane will usually allow the trajectories to describe successive 
states that the system may take as time proceeds but will usually not 
indicate the times at which the various states exist. For example, 
trajectories of displacement against rate of displacement would be suf- 
ficient to describe the successive states of motion of an oscillating 
system having one degree of freedom. The graphical representation of 
trajectories in the phase plane is facilitated by the construction of a 
family of isoclines where each isocline is a locus of points where the 
trajectories have a given slope (reference 15, pp. 10 and 2k8). Inves- 
tigation of the singularities in the phase plane indicates the type of 
motions which result from given initial, conditions, such as steady 
hunting oscillations or aperiodic or oscillatory responses that are either 
divergent or convergent. References I6 and I7 are two examples of applies 
tions of the techniques of nonlinear mechanics to cases in which an 
automatic -control con5)onent has rendered a system essentially nonlinear. 

The example chosen for this discussion is an aircraft having a non- 
linear variation of pitching-moment coefficient with angle of attack 
under the influence of angle -of -attack stabilization, (fig. -2(a)). The 
following assumptions are made: (1) the other coefficients are con- 

stants as listed in the appendix, (2) the reference angle of attack ai 
equals zero, and (3.) the nonlinear pitching-moment variation has a 
positive slope for smal 1 positive and negative values of angle of attack. 
As a consequence of the second assumption, the following equation will 
apply: 
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For the representation of the nonlinear function Cm(ct), consider 
the Cm curve in figure- 1 ha-vlng a positive slope at small value of a. 
If the a intercepts are symmetrical about a = 0, a cubic representa- 
tion of this curve would take the form 

Cjjj(a) = Cga3 - (12) 


where Cq and cq^ are constants. 


Considering a simplified notation for the equations of motion and 
substituting equations (11) and. (12) result in 


O O 

aj^D 6 + a^9 - c^a^ + C2^a = a^( -K]_a) 
b]^D9 + b2Da + b^o. = 0 


( 13 ) 


where Cl^ and Cj^^ have been assumed to be zero. The complete defi- 
nition of constants is given in the appendix. The equations are re-writ- 
ten as follows to illustrate that they can be reduced to an independent 
system: 


dO -^2 X 3 Cl 

— = 9 + — a-^ + 

d u a^_ V a^^ 


ai/ 


(I4a) 


da • 

dt " ^ ^ 


a 


(Ub) 


where ® ^ use. is made of the equation b2 = -bi. These equa- 

tions suggest that the complete motions can be examined by the instan- 
taneous states of a and 9. This idea is also apparent from the 
Laplace transform procedure suggested in a previous section since 
initial values of a and 9 are ^1 that are required for a definition 
•of the motion. Therefore, the a,0 plane is the phase plane for these 
considerations. 

Before actual trajectories are dra-wn in this plane, a family of 
isoclines is determined; The slope of any path in the phase plane 
is d©/da and is' determined by the division of equation (l4a) by 
equation (ikb). The resulting equation is 
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^0 3 

+ — a + 

(--■ 



ai 

1 

^1 / 


da 


-b3 • 

e + a 
b2 


(15) 


If d©/da is assigned a constant value X, equation (15) can be 
reduced to the form 


9 = 


£o ,3 / M X - Si - ^ a 
fO \b2 ai ai j 

^2 

\ + ““ 


(16) 


which is the isocline representing the locus of all points in the phase 
plane where trajectories must have the slope X. 


The singularities of the system axe defined as the equilibrium 
positions (steady-state values) that the system may attain arid are 
determined from equations (l4) by using 


d0 _ ^ _ 
dt " dt " 


(17) 


When the numerical values for the constants as given in the 
appendix are used, the following equation results for the isoclines; 

. - 309 (lO) 3 a 3 + ( 4 . 51 X + 257)a 

® X~4-"t.02, 


where X is the slope in radians per second per radian. The singulari- 
ties are found from equations (l4) when conditions of equation (I7) are 
imposed and are 

(1) a = 0, e = 0 

(2) a = 0.0278 radian,' 9 = 0.12^k radians per second 

(3) a = -0.0278, 9 = -O.I25U 

Figure 4 is a plot of the a,0 phase plane and shows some isoclines 
for various X. It is important to. note that the trajectories or paths 
which describe the motion have a definite direction that depends on the 
position in the phase plane. This directiori is determined from equa- 
tions (l4) by considering the value of equation (l4a) to be a vector 
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component in the 0 direction and the value of equation (l4b) to be a 
vector component in- the a direction. The direction of the resultant 
of these two components is the direction of the trajectory at that point 
in the phase plane j this resultant vector is termed the phase velocity 
at that point since the vector components determined from equations (l4) 
are themselves velocities. The arrows shown on the isoclines are not 
phase velocity vectors but indicate only the direction of the phase 
velocity on the various isoclines. 

Several trajectories are shown in figure 4 for .various initial con- 
ditions; these were drawn by the graphical method discussed in refer- 
ence 15 pages 248r252. It is evident from these trajectories and from 
the directions of- phase velocity on the given isoclines that the two 
singularities (±1.59°> ±0.1254) are points of stable equilibriijmj that 
is, the motion following a disturbance will oscillate about and tend to 
one of these stable states. A singularity of this type is teimed (refer- 
ences 15 and l6) a "stable focus" which means that the motion tends to 
this equilibrium in an oscillatory manner. The other singularity at 
(0,0) is an unstable equilibrium point. This instability is evidenced 
by the fact that the slightest disturbance from this point results in a 
motion to one of the stable foci. This singularity at the origin is 
termed a "saddle point." 

In this brief discussion it is evident that much can be learned 
about the characteristics of the system without resorting to actual 
solution of the nonlinear differential equations. The results presented 
for this case from an inspection of the nature of the phase plane can 
also be verified by the application of the various theorems presented 
in references l4 and 15 , which allow a determination of the stability 
of singularities and of such an important property as the existence or 
nonexistence of steady oscillations in the system. The methods of non- 
linear mechanics, however, are relatively new and, in their present state 
of development, they cannot be expected to handle an extreme variety of 
problems. The methods become quite involved even for systems having 
two degrees of freedom in which definition of the motion requires more 
than two quantities; the methods become almost impossible when more 
degrees of freedom are introduced. 


ANALOG RESULTS AND DISCUSSION 


As a part of the investigation of the effects of nonlinear stabi- 
lity derivatives on the transient longitudinal motions of an aircraft, 
the problem was studied through the use of electrical analog equipment. 
The analogs employed were the Engine and Control Simulator manufactured 
by Philbrick Researches, Inc., and the Reeves Electronic Analog Computer 
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(hereineifter designated the Philbrick Analog and the REAC, respectively), 
the units being, respectively, at the Lewis Flight Propulsion Laboratory 
and the Ames Aeronautical Laboratory of the NACA. Descriptions of the 
analog units are available from various sources, for example, refer- 
ences l8 and 19. 

The nonlinearities considered in these studies included both non- 
linear functions represented in equations (4), that is, the nonlinear 
variations of pitching -moment coefficient and lift coefficient with 
angle of attack. In all problems presented for both analogs, the non- 
linear functions were represented by straight-line segments although the 
methods of representation used in the two analogs differed. The other 
basic constants used are tabulated in the appendix, and the nonlinear 
functions are defined when presented. The presentation and discusslon- 
of results are given in two sections: (l) aircraft transient responses 

and (2) closed-loop transient responses. The first section contains 
analog solutions of the transients in angle of attack (and, in some 
cases, rate of pitching) resulting from step deflections of the canard 
elevator control surfaces. The second section presents some solutions 
for the closed -loop cases shown in figure 2. In both sections various 
types of nonlinear functions were chosen as being possibly representa- 
tive of the nonlinearities likely for the assiimed canard configuration. 

The results from the Philbrick Analog were obtained by photographing 
an oscilloscope; the REAC results were made with recording elements. The 
intent of the Philbrick studies was to survey qualitatively a wide range 
of nonlinear functions by observing and recording the more unusual tran- 
sient conditions. For this reason high accuracy in aircraft simulation ' 
was not required and the results are discussed in this light. In the 
REAC problems more exact solutions were required, the accuracy in 
final a and 9 transients being estimated at approximately 2 percent 
and the accuracy of individual REAC components being considerably better. 


Aircraft Transient Responses 


The analog results presented herein are aircraft transient 
responses in angle of attack (and pitching velocity) due to step deflec- 
tions of the canard elevator surfaces. Although only a few nonlinear 
situations are discussed in detail, the remarks which sxjmmarize the 
paper are based on a much larger number of cases for which analog 
solutions were obtained. 


Case I. Nonlinear Cjn(a): 


da ■ 


a < -2°, 


dCjn 

^ = nkj^, 0 < n < 1, 


-2°< a< 2°: 


a > 2°,- In case I the nonlinear function is the 
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pitching-moment coefficient due to angle of attack where the. slope at 
values of a less than- -.2° and greater than 2° is larger than the 
value of the slope in the region -2° < a. < 2°. In this case, the slope 
is always negative and stability is assured. Since the frequency of 
the transient oscillations is determined principally by this slope, 
changes in the frequency (period) are expected as the region in which 
the value of the steady- state angle of attack changes because of various 
elevator deflections. Little change in the damping envelope is 
expected. 

Results from the Philbrick Analog are presented in figure 5, 

In this run the value of kq was equal to "the slope 

for -2° < a < 2° is 0.463kq = -2.38. The a transient shown in 
figure 5(a) is the completely linear case where Cin(a) = kqa. ■ The total 
transient time shown represents 1.325 seconds. When the described non- 
linearity in Cia(a) is introduced, the results are shown by the four 
remaining records of figure 5. Curves of figures 5(b.) and 5(d) are 
two a transients, whereas curves of figures 5(c) and 5(e) are the 
corresponding nonlinearities. These nonlinearities were photographed 
by using the output of the element generating the nonlinear function as 
the vertical coordinate of the oscilloscope and the angle of attack as 
the horizontal coordinate. Hence, the region of angle of attack covered 
in the transients of figures 5(b) and 5(d) is shown again by the region 
of the nonlinearity pictured. The large spots shown on the nonlineaxi- 
ties represent the trim values of a due to the old and new 6 input. 
Hence, figure 5(c) shows that the values of & change in such a way as 
to make the trim values of a change from one extreme of the nonline - 
axity to the other. Since the slopes at these outer values of a are 
the same as that shown for the linear, case, the period should remain 
practically unchanged (compare figs. 5(a) and 5(b)). 

Figure 5(d) is the a transient when the change in 5 causes the 
trim a to change from a point where the slope is large to one in the 
region of reduced slope. The ratio increase in period (reduction in 
frequency) shown in figure 5(d) is approximately equal to the theoretical 
variation of being inversely proportional to the square root of the 
decrease in slope. (This variation will apply for lightly damped 
systems.) For the three transients shown, the damping envelope is seen 
to be essentially unchanged, the time required to reach steady state 
being about 1 second in each case. 

Two solutions from the REAC for this case (reproductions of actual 
recordings) are presented in figure 6. For these solutions the value 
of k]^ was -6.0 and n = 0.5^ and the cases shown are the results of 
moving the elevator in a programmed manner from 0° to about 8° and alter 
steady state is reached to about -k . The variations of 9', a, and 5 
with.'time are shown in this figure. The slight, phase difference between 
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control motion and the transients in the REAC solutions are due to the 
recording equipment and not to aircraft dynamics. The ±2° (±0.035 radian) 
positions noting the break points in the nonlinearity axe also Indicated 
on the a transient. The same general characteristics are observed 
for these two transients as mentioned previously, namely, the change in 
period and approximately the same time required for both transients to 
reach steady state. The steady-state values for a and 0 check rea- 
sonably well with those predicted theoretically. In figure 6 the 
responses show clearly that when the angle -of -attack variations are 
such as to be in the region of breaJts in the nonlinearity, the perio- 
dicity is not constant} when they, however, are in ein a region where 
the slope of the nonllneaxity is constant (near steady .state), the 
period does remain constant. For cases such as this illustration where 
it is relatively certain that the system is very lightly damped, this 
change in periodicity and finally a constant period near steady state 
can be attributed to a nonlinearity in the effective spring constant of 
the system, in this case Cni(a). 


dCm 

dq 


Case II, Nonlinear Cni(a): = ki, a < -2°, a >2°; 

= -nkj^, 0 < n < 1, -2° < a < 2°.- In case II the variation of the 


pitching-moment coefficient with angle of attack is nonlinear and con- 
sidered to have a positive (unstable) slope in the range - 2 ° < a < 2 °. 
A trim angle of attack for a given elevator deflection ( 6 ) cannot be 
expected to exist in the region - 2 ° < a < 2 ° since in this angle -of - 
attack region the transients are divergent. (This condition is shown 
by the presence of two real roots of opposite signs of the quadratic 
characteristic equation of equations (!»■), the negative root indicating 
a damped subsidence and the positive root indicating the divergent 
aperiodic mode.) 


Figure 7 presents some Philbrlck solutions for this case where 
ki = - 0 . 838 , n ^ 1. The total transient time shown represents 2.13 sec- 
onds. The a transient shown in figure 7(a) is for the completely 
linear case with CjQ(a) = kqa. The period of the damped oscillation 
is 0.284 second. Figures 7(l3), 7(d), and 7(f) are a transients for 
different 5 conditions, and figures 7 (c), 7 (e), and 7 (g) are the parts 
of the nonlinear function involved in each transient. 


From the nonlinearity shown in figure 7(c) the change in trim angle 
of attack' is large with the new trim angle of attack well into the 
region of stable slope. Therefore, little change is expected or noted 
in the transient of figure ' 7 (b) as compared with the linear case in 
figure 7(a-)- As the value of 5 is changed and the new trim angle 
of attack comes nearer to the break in the nonlinearity (see fig. 7 (e)), 
the initial portion of the transient of figure 7 (d) has become distorted. 



18 


NACA EM L50L11 


The characteristics of the oscillations about the new trim angle of 
attack, however, remain unchanged. As the 6 is changed only sli^tly, 
a new trim angle of attack is found to the left of the region of 
unstable slope as shown in figure T(g). The first overshoot of this 
transient (fig. 7(f)) is to an angle of attack in the unstable region, 
but the energy of the system is such that the complete unstable region 
is not traversed and the oscillatory characteristics about the new trim 
position become established. 


Only a very sli^t change in the 6 was found to be required for 
the change in trim angle of attack exhibited in the transients of 
figures 7(d) and 7(f) • Only when the new trim angle of attack was near 
the breaks in the nonlinearities did the transients indicate any peculi- 
arities which mi^t give an indication of the presence of a nonlinearity 
in the system. 

The REAC solutions presented in figure 8 for this case are for a 
value of k]_ = -3-0, n = 0.5» The two examples shown start from the 
same values of 6, trim a, and 9 with the firgt changing from 6 4° 

tO' S -2° and the second changing to 6 ~ -4 . In both of these 
conditions, the angle -of -attack variations are from a point of stable 
equilibrixom at positive a to a stable trim position at negative a. 

The transients indicate that the initial portions give the only hint of 
nonlinearity in the system. As noted previously, this initial part of 
the a transients becomes more out of the ordinary as the new position 
of trim a approaches the break point of the nonlinearity. These 
slight slope changes in the a transient are almost unnoticeable in 
the second curve. Examination of the 9 transients, however, indicate 
clearly the effect of the nonlinearity on the rate of pitching in the 
initial phase of the response. After this first effect on the tran- 
sients, the oscillatory characteristics become very regular. 


Case III. Nonlinear C-r(a): 


dC]- 

da 


= k, 


2 ^ 


a < -2" 


a > 2°j 


dCL 

5S- ' 


0 < n < 1, 


-2^ < a < 2^.- In case III the variation of 


lift coefficient with angle of attack is the nonlinear function having 
a larger slope at the outer values of a. Since the aircraft is lightly 
damped^ the primary contribution of the lift-curve slope to the oscil- 
latory characteristics is in the damping of the motion, Since^ for this 
case^ the slope of the pitching-moment -coefficient curve against angle 

/dCm \ ' 

of attack is a constant = ^1/^ change is to be expected in the 

period of the transient oscillations. ^ 


In figure 9 three Philbrick solutions are shown. Figure 9(a) is 
the completely linear case with k 2 = 3*^9 (li = ij and k]_ = - 1 . 90 - 
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The transient time is 2.48 seconds, the period is 0.194 second, and the 
cycles to damp to one-tenth amplitude are 3.39. Figure 9(h) is the 
transient which results when the trim a changes from a region of large 
slope (dCh/da) to a region of a of reduced slope. The value of n 
is 0.328 in this instance, that is, the slope at central a is about 
one-third of the slope at outer a. The period is essentially unchanged 
although the number of cycles to damp to one -tenth amplitude has 
increased to 4.20. For a change in trim a from the region of reduced 
slope to the region where the slope is the same as the linear case, the 
transient in figure 9(c). is practically the same as the one shown in 
figure 9(a) for the linear case. 

The change in slope used in this example is much greater than would 
be expected in a practical casej thus this change indicates that the 
actual effects of this parameter might result in very little change in 
the a transient characteristics. It is also emphasized that a simple 
examination of the transient (or even a comparison between figs. 9(b) 
and 9(c)) could not indicate that the nonlinearity was in lift-curve 
slope because the effect of variation of damping in pitch Cmq would 
also be manifested by changing the damping but not affecting the period. 
In actual flight tests, however, normal- acceleration- records might help 
to isolate this effect as a nonlinear Cl(oc). 

The REAC solutions sho-wn in figure 10 have k 2 = 4.19 and 
n = 0 . 834 . The linear equation of Cm(a) = -3.0a was also used in these 
runs. The variations of 5 shown were from '5=0° to 6 ~ 8° 
to b ^ -k . The periodicity is constant for these cases even though 
the variations of a in the second run are about the break point 
(a = -O.O 35 radian) of the nonlinearity. Since the two slopes of the 
nonlinear CL(a) function are nearly the same, the damping is practically 
unaffected in the two cases; in fact, plots of the logarithm of the peak 
amplitudes about steady-state values against time for these two cases 
indicate essentially the same slope and, therefore, the same damping, 
although the second run has a slight amount of scatter in the points. 
These cases which have a more reasonable value for n than used in the 
Philbrick solutions indicate the slight effects of this parameter under 
the conditions investigated. 

Case IV. Nonlinear Cm(a) and CL(a) .- In case IV the two non- 
linearities used in the REAC solutions for cases I and III were used 
together in the same run for two transient solutions. .These results 
are shown in figure 11 for 6 changes of 0° to 4° and 0° to 8°. In 
the first transient only the first overshoot is greater than a = 2°, 
the break point of the nonlinearities; hence, the transient character- 
istics are those defined entirely by the slopes of the nonlinearities 
at small a. In the second transient the period becomes constant after 
the first two cycles and exhibits characteristics similar to that of the 
first transient of figure 6 with the single nonlinear function Cm(a). 
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These results might have been expected since the nonlinear CL(a) had 
so little effect in case III. Other runs for cases where two nonline- 
arities were in the system indicated that the most noticeable changes 
in the transients occurred with changes in the Cm(a) function, as 
indicated previously. 

Review of effects of nonlinearities on aircraft transients . - Non- 
linear variations of the pitching moment resulted almost solely in 
effects on the period of the treinsient oscillations. Differences in the 
period of the oscillations in the different regions of angle of attack, 
due to the different slopes of the nonlinear function in these regions, 
were an effective way of noting the presence of this nonlinearity. 

Another noticeable effect was the irregularity in the period of the 
oscillations when the trim angle of attack for the various elevator 
deflections was near a break point of the nonlinearity. When the trim 
angle of attack is removed from the break in, the nonlinearity, the 
oscillation characteristics become very regular, that is, constant 
period and damping. If only moderate changes in trim angle of attack 
result from the control motion and the nonlinearity or a part thereof 
is traversed, the initial sections of the transients may indicate the 
presence of the nonlinearity through small slope variations, although 
the final oscillatory characteristics are quite regular. For large 
changes in trim angles of attack under the conditions noted, the small 
variations in the initial sections of the transient may become so 
slight as to be indistinguishable from the perfectly linear case. 

For cases when the nonlinear variation was in lift coefficient, 
large slope changes were required before changes in the oscillatory 
characteristics became appreciable for trim angles of attack in the 
different regions. This nonlinearity affected only the damping charac- 
teristics of the motions, the period remaining constant. This effect 
could be expected since very large changes in damping would be required 
to alter the period of this lightly damped system. 

When the two nonlinearities were considered simultaneously, with 
reasonable variations in each, the predominant changes in transient 
characteristics were similar to those noted for the single nonlinear 
pitching -moment variation, that is, the changes in periodicity previously 
mentioned with the time to damp remaining essentially unaltered. 


Closed-Loop Transient Responses 

The results in this section are transient response curves illus- 
trating the dynamics of the two closed-loop systems depicted in figure 2. 
The discussion will be in two parts: the closed loop with angle-of- 

attack feedback and the attitude stabilization, and includes Philbrick 
and REAC solutions for several cases. The nonlinear characteristics 
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will be those previously described for the aircraft) the error-sensing 
devices and proportional servomotors are considered as zero-lag systems. 
The transients are the responses to step input, command signals in ai 
and ©i . 

Angle -of -at tack feedback .- The block diagram for this case is shown 
in figure 2(a). Transients from the Philbrick Analog are shown for two 
cases: one where the nonlinear function is Cl ( a) and the other for the 

nonlinearity . C]n(a) . REIAC solutions were made for two different non- 
linear functions of pitching -moment coefficient. 

In the Philbrick solutions the value of Ki was approximately 
unity. For the case with nonlinear CL(a), the conditions are the same 
as described for case III of the previous section where the lift-curve 
slope at outer values of a is greater than the slope at small values 
of a. The transients shown in, figure 12 resulted from changes in 
that caused the final values of a to be in regions where the lift- 
curve slope is different. Figure' 12(a) shows the a transient about a 
steady-state value of a. in the region where the lift -curve slope is 
reduced. The oscillations In figure '12(b) are about a steady-state 
value of a in an a region of larger lift -curve slope. The increase 
in lift-curve slope is seen to improve the damping characteristics of 
the closed-loop system with angle-of -attack feedback in much the same 
way as the aircraft transients were improved. Investigation of the 
characteristic equation of this closed loop, the denominator of equa- 
tion (6), shows that the principal effect of such a change in lift- 
curve slope would be most evident in the damping. Comparison of tran- 
sients of figure 12 as well as an examination of the characteristic 
equation indicated that the frequency of the oscillations is constant. 

The Philbrick solutions with nonlinear Cni(a) were made for the 
case where the slope of the pitching-moment-coefficient curve was posi- 
tive (unstable) at small ct and stable at outer co. The outer slopes 
were -1.90, and the center slope was 0.655 . Two a transients illus- 
trating the effects of this nonlinearity on the closed-loop responses 
are given in figure I 3 . In figure 13(a) the change in reference angle 
results in a change in the steady- state values of a from one 
region of stable slope to -the other. In figure 13(b) the input change 
causes the final value of a to be a value in the region of angle of 
attack where the aircraft has an unstable pitching -moment curve, and 
the transient is satisfactory. The physics of this change is explained 
by the action of the angle -of -attack feedback causing the control sur- 
face deflection 6 to give a moment proportional to the instantaneous a 
thus, this moment adds to the effective spring constant of the system. 

In this instance, the value of Ki is large, enough so that this con- 
tribution overcomes the destabilizing effect of the spring in the air- 
craft at central a and renders the over-all system ■ stable in this 
region of a. This stability could also be determined from the 
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inequalities (7) which show that the value of Ki used assures com- 
plete stability of this closed-loop system. 


REAC solutions were also made for a case where the pitching-moment 
curve had an unstable slope at small values of angle of attack. The 
same conditions as discussed for REAC runs in case II were used here 


with a value of 


dCm 

da 


-3..0 for the stable slopes at outer a and 


dCm o o 

= 1.5 ' for -2 < a < 2 . The value of' Kq was 1.0. The three a 

transients of figure Ik are for the , step inputs of .1^, k^ , and 8^, 
The 2 position noting the break in the nonlinearity is shown on each 
curve. For this case^ an examination of the inequalities (7) reveals - 
that the value of Ki = 1.0 is not large enou^ to assure complete 
stability in the region -2° < a < 2°> hence, no steady- state value 
of a ^is possible in this region. Therefore, for even the small input 
ai = 1 , the final value of a is greater than the break point of the 
nonlinearity. The measured values of steady-state a due to the inputs 
involved are 2.49°, 3*21°^ and 4.24° con^pared to the theoretical values 
of 2.50 , 3-26^, and 4.30°> respectively. Since the value of a is 
greater than 2 after the first cycle for all the transients shown, the 
period of the oscillations is essentially the same for all cases. 


The final REAC solutions are shown for the case where the Cm(a) 

dC_ 

nonlinearity is the same as case Ij namely, = -3.0 where 

dC 

“2 < a < 2° and = -6.0 when a < - 2°, a > 2°. The inputs 

were again 1°, 4°, and 8°, and the value of Ki was 1.0. The resulting 
closed-loop a transients are given in figure I5.. In these runs too, 
the steady-state values of a check satisfactorily with the theoretical 
values. Note that figure 15(c) for ai = 8^ has a final value for a 
of about 2°; hence , the oscillations, are about the break point of the 
nonlinearity. The period of this transient is very constant and is 
different from the period of figures 15(a) and 15(b) which both remain 
in the region 0 < a < 2°. The period of figures 15(a) and 15(b) is 
about 1.17 times the period of figure 15(c). A theoretical value of 
the change in period, obtained from consideration of the values defined 
by the two slopes of Cm(a), is 1.15j this value is close to the meas- 
ured value even thou^ the oscillations were not in the region where 
this result could be expected to apply. Stability of the combination of 
the aircraft and autopilot with control proportional to angle of attack 
was obtained by adjustment of the autopilot proportionality constant, 
and the effects of the nonlinearities were similar to those noted for 
the aircraft transients to step control deflections. 
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Examination of figures l4 and 15 indicates that all the cases shown 
with nonlinearities in the system exhibit a steady-state error between 
the reference ai and the actual a. This characteristic is not due 
primarily to the nonlinearities of the system. An investigation of a 
linear aircraft in this system of angle -of -attack feedback reveals that 
the complete system is not a zero-position-error system; therefore, 
unless some form of compensation is included in the loop the output will 
not equal the input at steady state. For the cases shown, however, 
adequate compensation of the system with nonlinearities as described 
may be no more difficult than that required to make the linear system a 
zero-position-error-system. 

Attitude stabilization :- In the results (REAC) to be presented for 
the problem of attitude stabilization", the only nonlinearity considered 
was Cni(a). In the first example the Cm(a) function used for REAC solu- 
tions for case I, dCm/da being negative for all a but having a larger 
slope at outer a, was employed. The second example is the Cin((x) func- 
tion for case II, dCm/da having a stable (negative) slope .at outer a 
and an unstable slope at -2° < a < 2° and a K2 which does not assure 
stability at central a. The analysis in previous sections indicated 
that the first example should exhibit constant steady-state values after 
transients die out, whereas the second example should show continuous 
hunting oscillations. For both cases K2 = 1.0. 

Figure 16 illustrates the first example where the Cm(a) has a 
negative slope at all a. Figure l6(a) is for = 4.6° and fig- 
ure 16(b) is for 0i = The upper curves in figure l6(a) and fig- 

ure 16(b) are the a variations, whereas the lower curves are of the 
error (9i - 0). In both of these cases the error tends to zero indi- 
cating the zero position error characteristic of the system. The cubic 
characteristic equation of this system, see equation (9), would show an 
oscillatory mode of motion superimposed upon the convergent aperiodic 
mode and is indicated by the negative real root. From the characteris- 
tics of the system at central a,, the oscillatory mode damps to one -half 
amplitude in about one -third of the time required for the aperiodic mode 
to damp the same amount. In figure l6(b), the period becomes constant 
only after the variations of a become less than 2°, whereas in fig- 
ure 16(a) the period is constant throughout. 

The results for the second example where the Cm(a) function had an 
unstable slope in the region -2° < a, < 2° are presented in figure I7. 
Curves in figures 17(a), 17(b), and 17(c) were results of step inputs 
of ©i of 0.7°, 4.6°, and 8.7°, respectively. In each set the time 
variations of a and (0^ - 0) are shown. Though the initial portions 
of the responses are different for the three cases, the resulting 
hunting oscillations are essentially the same. These curves show the 
constant -amplitude constant -period oscillations predicted in the 
theoretical analysis for this type of nonlinearity. 
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Stability of the attitude -stabilization system, when the control 
surface deflection is proportional to the error in attitude,- could be 
completely assured for the types of nonlinearities considered herein 
with proper adjustment of the autopilot proportionality constant. The 
existence of continuous hunting oscillations for a particular nonlinear 
pitching-moment variation and improper autopilot adjustment was predicted 
theoretically and demonstrated by analog solutions. The pilot of a 
conventional aircraft might apply control in a manner that is approxi- 
mately proportional, though probably with a time lag, and such a hunting 
oscillation as described mi^t be realized if that aircraft had a charac- 
teristic like the nonlinearity in the data of figure I 7 . 


CONCLUSIONS 


This study has been concerned with the effects of two nonlinear 
stability derivatives, caused by the nonlinear variations of pitching- 
moment and lift coefficients with angle of attack, on the longitudinal 
motions of an aircraft. In addition to the theoretical methods dis- 
cussed, the excellent adaptation of analog equipment to this dynamical 
study gave results for the responses of a canard aircraft to step 
control deflections and to the influence of two types of proportional 
automatic control. From the theoretical considerations and the analog 
results presented herein, the following conclusions summarize the 
investigation: 

1. When the nonlinear functions are approximated by a series of 
linear segments, the methods of the Laplace transformation allow the 
transient responses of the aircraft alone and also with proportional 
automatic control to be calculated. The computational procedure 
required, however, is tedious and time consuming. 

2. The most noticeable effects of the nonlinear pitching -moment 
variations on the response of the aircraft to step control deflections 
was in the periodicity of the oscillations . The occurrence of an 
irregular period and of different constant -period oscillations about 

the various trim angles of attack was the most satisfactory way of noting 
the presence of this nonlinearity. 

3* Nonlinearity in the lift curve affected only the damping of the 
transient angle -of -attack oscillations, the period remaining constant. 
Large slope changes in this nonlinear function were required before the 
effect on the time to damp became appreciable. 

4. Simultaneous occurrence of the two types of nonlinearities con- 
sidered herein, with reasonable values of each, resulted in transient 
characteristics which were similar to those with the single nonlinear 
pitching-moment function. 
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5. StalDility of the combination of the aircraft and autopilot with 
control proportional to angle of attack was obtained by adjustment of the 
autopilot proportionality constant^ and the effects of the nonlinearities 
were similar to those noted for the aircraft transients to step control 
deflections. 

6 . Stability of the attitude-stabilization system^ when the control 
surface deflection is proportional to the error in attitude, could be 
completely assured for the types of nonlinearities considered herein 
with proper adjustment of the autopilot proportionality constant. The 
existence of continuous hunting oscillations for a particular nonlinear 
pitching -moment variation and improper autopilot adjustment was 
predicted theoretically and demonstrated by analog solutions. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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APPENDIX 

FACTORS AND COEFFICIEHTS USED IN TBE ANALYSIS 


The complete expressions for the factors and coefficients used in 
the analysis section of the paper, in the order in which they appear 
in the text are as follows: 





b = 


aibg 


®-l 


|l.C.(s^ 


a 


a(0)s + 


a2 + ai^ 
_ 


a(0) + 0(0) 



^1^2 


2Cc% = 


ag a4 k2 

+ + r— 

ai ai hi 



a2k2 

a^^b^ 


ai 


2 

(When ( 0 ^ >0 and 2^c% >0, a)n corresponds to the undamped natural 

frequency of the short -period mode of motion of the aircraft alone 
and ^ represents the conventional nondimensional damping ratio.) 


^5 

a^^ a^^bg 


d = 




+ ^iH 
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0 ( 0 ) 


[l.C.(s)]^ = e(o)s^ + 

0(O)^^j + a(0) I- 


^2 ^ 


— + — + >r + ®(o) 

ai ai hj 


s + 0(0) 


agk^ 


kg\ / a 


l^kg ki 




^2^ ” 

aibg 


+ + ki T) 

s\ ^1^2 


= :;sF 


qSc 


a2 = 


c S— 
■^“q 2V 


aij. = 


-C — 
2V 


^5 = 


>^ - K _ ®V 

- -^2 - ^ 


^4 = Cl 

o 

■ Values used for example of method of nonlinear mechanics are as 
follows ; 

Co = -546. 

Cl = -1.5 

ai = 0.00177 

a2 = 0.00712 

a4 = 0 

a^ = 1.045 

hi = -b2. = 0.774 
b^ = 0 

dCr 

b3 ^ ^ = -3.49 

J da 

Ki =1.0 . 
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For the solutions made on the Philbrick and J^eeves Analog equip- 
ment^ the numerical values for the aforementioned constants a^^, a2, 

a^, a^, bj^, and b|^ were used. The values for these coeffi- 

cients resulted from the following flight conditions and stability- 
derivatives: 


Mach number = 1.8 
V =2009 ft/sec 
q = 4800 lb/ft2 

s = 2.52 ft^ • 
c = 1.4 ft 
ly' = 30 slug -ft ^ 
m = 4.66 slugs 


^m ' 

= - 20.43 

2V, 


dc. 


'^m , 

^“6 = ^ 


"Lb = " 
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Figure 1.- Typical nonlinearities in pitching-moment and lift coefficients. 
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(a) Angle- of-at tack stabilization. 



(b) Attitude stabilization. 

Figure 2.- Block diagrams illustrating the two types of proportional 
automatic control under consideration. 
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Figure 3.- Illustration of the method of approximating nonlinear functions 

with linear segments. 
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(b) a transient. 


(c) Nonlinearity. 



(d) a transient. 

Figure 5*- Philbrick solutions for case I. 


(e) Nonlinearity. 

Aircraft a transients and 


parts of nonlinearity involved; nonlinear 
da 


dC 


= - 5 . 14 , 


a < -2°, 


a > 2°; 


= - 2 . 38 , - 2 ° < a < 2°; — = 3 .^ 9 , all a. 

da 


da 
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(b) a transient. 



(d) a transient. 



(f) a transient. 
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(c) Nonlinearity. 
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(e) Nonlinearity. 
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(g) Nonlinearity. 


Figure 7.- Philbrick solutions for case II. Aircraft 
parts of nonlinearity involved; nonlinear Cjn(a), 


dC 
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dCr 


a 



da 


transients and 

= -0.838, 
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1^3 



(c) n = 0.328. 


Figure 9.- Philbrick solutions for case III. Aircraft a transients; 


nonlinear 0^(0); ^ = 3.49, a < - 2 °, a > 2 °; ^ = n( 3 . 49 ), 


dCL 


dC., 


da 


da 


-2° < a < 2°; — £ = -1.90, all a. 
da 
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(b) Steady-state a in region 
dCr 

where — = 3«^9* 
da 


Figure 12.- Philhrick solutions for proportional control system with 


angle-of-attack feedback, a transients, nonlinear Cj^(a), = 3.^9 


a < -2“, a > 2°, ^ . i.u, -2“ < „ < 2<>; ^ . .1.50, all a, 


dC. 


da 


autopilot constant % 1.0. 
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Figure I3.- Philbrick solutions for proportional control system with 

angle-of-attack feedback, a transients and parts of the nonlinearity; 

nonlinear Cj„(a), ^ -I.90, a < - 2 °, a > 2 °; ^ = O.655, 

dCi 

-2 < a < 2 ; ='3.49^ all a; autopilot constant a: 1.0. 





PAGE MISSING FROM AVAILABLE VERSION 




52 


NACA RM L50L11 



OC, 

f*€idtan 



(b) tti = 4°. 



(c) = 8 . 


Figure 15.- REAC solutions for proportional control system with angle-of- 

dC 

attack feedback. Nonlinear Cjj(a), — S = - 3 . 0 , -2° < a < 2°: 

da 

^ = -6.0, a < -2°, a > 2°; ^ = 3 . 49 , all a; autopilot 
constant = 1.0. 



NACA RM L 50 L 11 


53 


oc 

radian 





(Gi-eX 

rad/on 








ssaskmaa 




(t>) e± = 9 ' 


oc, 

raJion 





(Gi-OX 

radian 




Time , , sec. 

(a) ©i = 4.6°. 

Figure l6.- REAC solutions for proportional control for attitude stabili- 
zation., Nonlinear Cm(a), ^ =-3.0, -2° < a < 2°; ^ = -6.0, 

o o dCi do, 

a < -2 , a > 2 ; = 3 . 49 ; autopilot constant K2 = 1.0. 


NACA RM L50L11 



3.^9j autopilot constant Ko = 1.0, 
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